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Abstract 

We analyze the 3-parameter family of exact, regular, static, spher- 
ically symmetric perfect fluid solutions of Einstein's equations (corre- 
sponding to a 2-parameter family of equations of state) due to Pant 
and Sah and "rediscovered" by Rosquist and by the present author. 
Except for the Buchdahl solutions which are contained as a limiting 
case, the fluids have finite radius and are physically realistic for suit- 
able parameter ranges. 

The equations of state can be characterized geometrically by the 
property that the 3-metric on the static slices, rescaled conformally 
with the fourth power of any linear function of the norm of the static 
Killing vector, has constant scalar curvature. This local property does 
not require spherical symmetry; in fact it simplifies the proof of spher- 
ical symmetry of asymptotically flat solutions which we recall here for 
the Pant-Sah equations of state. 

We also consider a model in Newtonian theory with analogous 
geometric and physical properties, together with a proof of spherical 
symmetry of the asymptotically flat solutions. 

1 Introduction 

Solutions of the Tolman-Oppenheimer-Volkoff equation yield quite realistic 
models for non-rotating stars. Moreover, the system of the static Einstein 
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equations with perfect fluids also provides a testcase for general mathematical 
techniques and has stimulated their development. Firstly, group theoretical 
and Hamiltonian methods for generating solutions have applications to this 
system on the local level (see, e.g. [HE])- Secondly, under natural global con- 
ditions like asymptotic flatness the system is overdetermined, which should 
lead to spherical symmetry of the solutions for all equations of state (EOS) 
p{p) with p > for p > 0. This long-standing conjecture has, in essence, been 
settled recently by Masood-ul-Alam p] by using extensions of the techniques 
of Witten's positive mass theorem [1]. Thirdly, as an interesting result on 
ODEs Rendall and Schmidt [5] and Schaudt [6] have proven existence of 1- 
parameter families of spherically symmetric, asymptotically flat solutions of 
the Tolman-Oppenheimer-Volkoff equation for very general classes of EOS, 
and for all values of the central pressure and the central density for which 
the EOS is defined. 

To illustrate these mathematical results, to study remaining conjectures 
(see e.g. [3 [8]) and to make the connection with physics, it is desirable to have 
at one's disposal some "exact" model solutions, preferably with physically 
realistic equation of state. Well suited for these purposes is the two parameter 
family of Pant-Sah EOS (PSEOS) [9] 

p = p_(l-A)^ + p+(l + A)^ (1) 
P = ^[P-(1-Ar-P+(1 + Ar] (2) 

for some constants A, p_ and p+ with < A < 1 and < p+ < p_. This is a 
parametric representation of solutions of the second order ODE I[p{p)] = 
where 

I[p{p)] = y + 2n+ip + p)^ with n=P±^^. (3) 
5 dp p + 3pdp 

Putting p_|_ = and eliminating A in ([1]), ([2]) we obtain the one-parameter 
family of Buchdahl equations of state (BEOS) with has the 2-parameter 
family of Buchdahl solutions [TU]. The case p = const, is included here as 
a (degenerate) solution of / = 0, and it also arises in the limit p+ — ^ p_, 
A in ([I]),® (c.f. Sect. 2.2.2). The general 2-parameter family (^,(^ 
was considered by Pant and Sah [9] who gave the 3-parameter family of 
corresponding solutions in terms of elementary functions. A decade later, 
the PSEOS arose in the course of work on a uniqueness proof [HI [12], which 
lead to the first "rediscovery" of the Pant-Sah solutions [13]. Moreover, the 
Pant-Sah solutions also came up in a systematic Hamiltonian approach to 
relativistic perfect fluids by Rosquist |2j|- 
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These papers also established the basic properties of the solutions relevant 
for their use as stellar models, namely: 

• The Pant-Sah solutions are regular as long as the central pressure stays 
bounded (contrary to the claim by Delgaty and Lake [I4j; see however 

for a correction). 

• All Pant-Sah solutions except for the Buchdahl solutions have a fluid 
region of finite extent, which is obvious from ([T]), since p > at 
p = iff p+ > 0. 

• The energy density is positive and the pressure is non-negative every- 
where, and these functions decrease monotonically with the radius. 

• Under suitable restrictions on the parameters, the speed of sound re- 
mains sub luminal everywhere |16j . 

Morover, Pant and Sah showed that the parameters can be fitted quite well 
to neutron star data [D], while Rosquist JW\ considered Pant-Sah solutions 
as possible "traps" for gravitational waves. 

The purpose of this paper is to give a unified description of the Pant- 
Sah solutions, to explain their wide range of applicability and to extend 
it even further. Apart from the physically relevant properties mentioned 
above, we find here the following. Depending on the choice of p+ and p_, 
the mass-radius relation (which is a polynomial equation quadratic in the 
mass) is either monotonic, or it exhibits a maximum of the radius only, or 
a maximum of the mass as well before it reaches a solution with a singular 
center. However, the surface redshift and therefore the quotient mass/radius 
uniquely characterizes a Pant-Sah solution for any given EOS. This implies 
that the mass-radius curves can form a single, open "loop" but will not 
exhibit the "spiral" form typical for degenerate matter at extreme densities 
[T71 [18] . Nevertheless, for suitable p+ and p_ the mass-radius curve fits 
remarkably well with some quark star models discussed in the last years 
(see, e.g. [l9]-[22]) except at extreme densities. 

As to the mathematical properties of the Pant-Sah solutions, the key 
for their understanding is the " Killing- Yamabe property". By this term, 
motivated by the Yamabe problem [23] , we mean the following: For all static 
solutions of Einstein's equations with a perfect fluid (defined only locally and 
not necessarily spherically symmetric) we require that gf- = (1 -|- fVYgij/lQ 
is a metric of constant scalar curvature TZ^ where gij is the induced metric 
on the static slices, V is the norm of the Killing vector and / is a constant 
chosen a priori [13]. (If / 7^ 0, it may be absorbed in V by a suitable 
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scaling of the Killing vector. The case / = clearly corresponds to fluids 
with p = const.). If the Killing- Yamabe property holds, the field equations 
imply that g^j = (1 — fVYgij /IQ, has constant curvature 71- as well, with 
TZ- 7^ TZ+ in general. Together with dH),® and / = 0, the Killing- Yamabe 
property provides a third alternative characterization of the PSEOS, and the 
two curvatures are related to the constants in ([T]) and ^ by 1Z± = 5127ip±. 

To understand how the Killing- Yamabe property leads to the "exactness" 
of the Pant-Sah solutions, we note that spherically symmetric 3-metrics with 
constant scalar curvature and regular centre are "Einstein spaces" (i.e. the 
Ricci tensor is pure trace). Such spaces enjoy simple expressions in suitable 
coordinates, and the same applies to the conformal factors defined above. 

To sketch the proofs of spherical symmetry of asymptotically fiat solu- 
tions, we need two further generally defined rescalings of the spatial met- 
ric, namely g*j = K(y)gij where the function K(y) is chosen such that g*j 
is fiat if gij is a Pant-Sah solution {K{V) is non-unique in general), and 
g'.. = (1 - pV'^)'^gij/lQV\ We then show that for solutions with PSEOS 
there holds the "±"- pair of equations 

A'7^, = P±Bpt (4) 

where A' is the Laplacian of g'^j, Bf^ and B~j are the trace- free parts of the 
Ricci-tensors of gf^ and g~j respectively, and /5+ and /?_ are non-positive 
functions. 

To show uniqueness of the Buchdahl solutions for the BEOS = 0, 
TZ- 7^ 0), one first shows that all asymptotically fiat solutions must extend 
to infinity [3, [211 ES] ■ Hence for a Killing vector normalized such that V ^ 1 
at infinity, we can choose / = 1 and K{V) = (1 + 1^)^/16 in the definitions 
above. Then there are two alternative ways to continue [H]. The first one 
consists of noting that g*j = gj^ is asymptotically fiat with vanishing mass. 
Hence the positive mass theorem [H [26] implies that these metrics are fiat 
and {v., gij) is a Buchdahl solution. Alternatively, we can integrate the "mi- 
nus" version of (jl]) over the static slice. By the divergence theorem and by 
asymptotic fiatness, B~j = 0, i.e. g~j is an Einstein space, which again yields 
a Buchdahl solution. 

In the generic case 71+ ^ 0, the divergence theorem alone applied to (jlj) 
is insuffient for a proof as TZ^, cannot be made on the fiuid boundary in 
general. However, by employing a suitable elliptic identity in the vacuum 
region as well, the maximum principle now yields 7Z^ > 0. Then the positive 
mass theorem leads to the required conclusion. 

The positive mass theorem combined with Equ. (HI) has been been em- 
ployed earlier in proofs of spherical symmetry in the cases of fiuids with 
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constant density [27j and "near constant density" ^28j. Moreover, there are 
generalizations of (jl]) for non-Killing- Yamabe EOS, for fluids which satisfy 
J < 0, which again give uniqueness fl2[ [29] . The general proof of spherical 
symmetry by Masood-ul-Alam [3] involves modified Witten spinor identities 
and integral versions of generalizations of (jll). 

Jointly with the PSEOS we will consider here a model in Newtonian the- 
ory characterized by the following counterpart of the Killing- Yamabe prop- 
erty: We require that a conformal resettling of flttt spaee with {v — v^jVo, 
where v is the Newtonian potential (not necessarily spherically symmetric) 
and V a constant, is a metrie of eonstant settlar eurvttture TZ-. This leads to 
the 2-parameter family of equations of state 



where IZ^ = 5127rp_ and is another constant which has here no obvious 
relation to curvature. We will refer to as "the Newtonian equation(s) 
of state" (NEOS). For p+ = the NEOS are polytropes of index 5 which 
are analogous to the BEOS. The general NEOS and the corresponding so- 
lutions may be considered as "Newtonian limits" of the PSEOS and the 
Pant-Sah solutions, with similar properties for low density and pressure. As 
to uniqueness proofs for asymptotically flat solutions with the NEOS, there 
is available some sort of counterpart of (jl]), and the positive mass theorem 
has to be substituted here by the "virial theorem". 

This paper is organized as follows. In Sect. 2 we give the field equations 
in the Newtonian and in the relativistic case and introduce our models. In 
Sect. 3 we rederive the spherically symmetric solutions and discuss their main 
properties, in particular the mass-radius curves. The Section 4 we prove 
spherical symmetry of asymptotically fiat solutions with the NEOS and the 
PSEOS. The Appendix contains general material on conformal rescalings of 
metrics and on spaces of constant curvature. 



Our description of the Newtonian and the relativistic fluids will be as close as 
possible. For simplicity we will use identical symbols {gij, g~j, Vj, TZ, TZ~ ...) 
for analogous quantities but with different formal definitions, depending on 
the context. 

We denote by JF the fluid region, which we assume to be open and con- 
nected, and which may extend to infinity. V is the open vacuum region 
(which may be empty) and dT = dV is the eommon boundary (i.e. infinity 




(5) 



2 The Field Equations 
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is not included in dV ). This redundant terminology is useful to describe the 
matching. When J-' is spherically symmetric it is called "star". 

2.1 Newtonian Fluids 
2.1.1 General properties 

We consider as Newtonian model JF U dJ^ U V a manifold {Ai,gij) with a flat 
metric Qij. The potential function v is assumed to be smooth in JF and V, 
C^'^ at (9JF, negative everywhere and f — at infinity. For smooth density 
p(x*) and pressure functions in JF, with p — at dJ^, Newton's and 

Euler's equations read 

Av = A-Kp (6) 
Vip = -pViV (7) 

where Vj and A = VjV* denote the gradient and the Laplacian of fiat space. 
(Indices are moved with gij and its inverse g^^). 

A general EOS is of the form H{p, p) = 0. (If possible we choose H{p, p) = 
p—p{p)). H{p,p) should be defined in the intervals p G [ps, oo) andp G [0, oo) 
with Ps > and smooth in the intervals p G {ps, oo) and p G (0, oo). Using 
the EOS, we can write p and p as smooth functions p{v) and p{v) and Euler's 
equation ([7j) as dp/dv = —p. 

To recall the matching conditions, we note that, from the above re- 
quirements, the metric induced on dJ-" is C^'^ and the mean curvature of 
dJ-" is continuous. We now write these conditions in terms of the quantity 
w = VjfV*f . (The gradient always acts only on the subsequent argument, 
i.e. w = (Vif)(V*f)). From ([HD it follows that the quantity in brackets on 
the l.h. side of 

[w-^V'vViW - 87rp] = [w^^V'vViw] (8) 

is continuous at the surface. Hence dH]) must hold, where and 
" dV <^=" denote the approach to the boundary from the fluid and the vacuum 
sides, respectively. 

Generalizations to several disconnected "matching surfaces" are trival and 
will not be considered here. 

To formulate the asymptotic properties we consider an "end" = 
A4\{a compact set}. We assume that, for some e > 

(9) 
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where M is the mass. With ^ and ([7|) this imphes that 

P = 0{^1 P = 0{^). (10) 

A more natural but involved precedure is to derive the falloff conditions of 
the potential v and of p only from the falloff of p (c.f. 0). 



2.1.2 The Newtonian equation of state 

To introduce our model it is useful to rescale the Euclidean metric by the 
fourth power of a linear function of the Newtonian potential, i.e. we define 
Qij = Qijiy " vY I^Q for some constant > 0. 

We use the general formula (|89ll for conformal rescalings with pij flat, 
$ = {v — v)/2, and hence pij = g~y Together with the field equation ([6]), we 
find that 

IZ^iy -vf = hl2Tip (11) 

where TZ- is the scalar curvature of g~y We now determine the NEOS by 
requiring that TZ- = const.. Introducing p_ by TZ- = 5127rp_ and another 
constant p+, eqs. ((Tj) and (ITT]) yield 

p = p-{v-vf p = hp_(^v - vf - p+]. (12) 

o 

We require that p+ > and p_ > 0. Eliminating the potential we obtain 
the NEOS equ. ([5]). In terms of the variables p/p- and p/p_ and r = 

{p+/P-Y^^ this equation reads p/p- = | (p/P-)^ — ■ This means that 
we have singled out p„ as a "scaling" parameter while r G [0, oo] plays a more 
"essential" role. (This terminology mainly serves to simplify the analysis of 
static spherically symmetric solutions in Sect. 3. Both parameters have 
direct physical significance, as follows from ([T3j) below. On the other hand, 
from the dynamical system point of view, both parameters can be considered 
as "scaling" except in the case p = const, (c.f. [HUl IHT] ). 

Fig. (pP) shows the NEOS for the values r = (3 - v^)/2 ^ 0.382, r = 0.6 
and r = [(-v/S — l)/2]^/^ 0.786. (These particular values play a role in 
Relativity and are chosen here for comparison). 

The speed of sound C defined by the first equation of (|T3l) takes the simple 
form 

in terms of the potential. We recall that v was taken to be positive, hence 
dp/ dp > and (fT3l) makes sense. 



7 





0.2 0.4 0.6 0.8 1 1.2 

P/P- 

Figure 1: The equation of state ([5]) for the values r = 0.382 (thin hne), 
r = 0.6 (medium) and r = 0.786 (thick hne). 

At the surface where the pressure is zero the potential, the density and 
the speed of sound take the values 

Vs = V-T, ps = T'^P+, Cl = ^T. (14) 

5 

Note that ps and Cg are determined by the equation of state alone, as 
opposed to Vg and v which will be used in Sect. 3.1 to parametrize the 
spherically symmetric solutions. The polytrope of index 5, for which the 
solutions extend to infinity, arises from the equations above as the special 
case p+ = Ps = = Vg = V. The corresponding curve would pass through 
the origin in Fig. ([T]), very close to the curve for r = 0.382. 

2.2 Relativistic Fluids 
2.2.1 General properties 

We consider static spacetimes of the from ]RxA^ = Mx(jFU dJ-" U V) with 
metric 

ds'^ = -V'^dt^ + Qijdx'dx^ (15) 

where V{x^) and gij{x^) are smooth on and V and C^'^ at dJ-'. Moreover, 
< \^ < 1 on and ^ — > 1 at infinity. On JF we consider smooth density 
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and pressure functions p(x*), p(a;*), with p ^ on (9jF, in terms of which 
Einstein's and Euler's equations read 

= 47rV(p + 3p) (16) 
7^., = \/-V,V,y + 47r(p-p)(7,, (17) 
Vip = -V-\p + p)ViV. (18) 

The gradient Vj, the Laplacian A = VjV* and the Ricci tensor IZij now refer 
to Qij. As well known the Euler equation f|T8l) is a consequence of the Bianchi 
identity for TZij. 

A general equation of state H{p,p) = 0, in particular with H{p,p) = 
P " p{p): should be defined in p G [ps, oo) and p G [0, oo) with ps > and 
smooth in the intervals p G (ps,oo) and p G (0,oo). Euler's equation ( |T8l) 
together with the equation of state imply that there are smooth functions 
p{V) and p(y), and Euler's equation becomes dp/dV = —V~^{p + p). 

In analogy with the Newtonian case, the metric induced on dJ-" is C^'^ 
and the mean curvature of dJ-" is continuous. In terms of the quantity W = 
ViW^V, the matching conditions together with equ. (fT6ll imply that the 
quantity in brackets on the l.h. side of 

[W-'V'VV,W - SnVp] = [W-'V'VV,W] (19) 

is continuous at dJ-' and hence (fT^ holds. 

To formulate the asymptotic properties we consider an "end" = 
A4\{a compact set}. We require that, for some e > 

Ml Ml 

d,d,V = -dAy + 0{^J, (20) 

2M 1 2M 1 

gij = (IH )Sij + 0{——), dkQij = dk 6ij + 0{ 



dkdiQij = dkdi—6ij + 0{^—) (21) 



fy ' i^3~|~f ' 



in suitable coordinates, where M is the mass. Equs. (1161) and f|T8|) together 
with the decay conditions fpUj) and fl?Tl) imply that 

P = 0(;:^)' P = 0(;i7)- (22) 

Here the falloff conditions of the potential V could also be derived from some 
weak falloff conditions of g^, and p and p [7]. Clearly a substantial refinement 
of all asymptotic conditions is possible if is vacuum, c.f. 
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2.2.2 The Pant-Sah equation of state 

We now introduce conformal rescalings of tlie spatial metric of the form 
9ij = 9iji^ ^ /^)^/16 for some constant / which we take to be non-negative 
(this just fixes the notation), and we restrict ourselves to the range V < 1/ f. 
While any / > could be absorbed into V by rescaling the static Killing 
field, we have already fixed the scaling above by requiring V" ^ 1 at infinity, 
which is why the extra constant / will persist here in general. We now use 
the standard formula (|89l) with p = gij, $ = (1 ± fV)/2 so that pij = gfy 
Together with the field equations (|T6l) and (|T7I) this gives 

^7^±(l ± fVf = -(A - ^n){l ± fV) = 27r[p(l T fV) T QfpV] (23) 

where TZ and Tl± are the scalar curvatures of g^j and g^, respectively. By 
differentiating (PHI) with respect to V we obtain 



dn± _ 25607r(p + 3p) 

~W ~ (1 ±/y)6 



where k has been introduced in ([H]). 

We now implement the "Killing Yamabe property" defined in the in- 
troduction by requiring that at least one of the curvatures TZ+ and TZ- is 
constant. This implies that the quantity in brackets in ( 12^ vanishes, i.e.. 



''1'^' (25) 



(1 - PV^) 

and therefore the other scalar curvature is necessarily constant as well. Using 
this in fl23p and setting TZ± = 512ttp± we obtain 

p = p4l-fV)' + p+{l + fVf, (26) 

P = g^[P-(l-/^r-P+(l + /^n (27) 

which yields the parametric form ([I]), (E]) of the PSEOS when we set fV = A. 
Positivity of the pressure now requires that we restrict ourselves to < p+ < 
p_ < oo. Note that Equ. ( |25l) implies that dp/ dp takes on all real positive 
values for the allowed range < fV < 1 of V^. 

The case p = const, is included in f l2B|) and flTH) in the limit p+ — > p^ 
and fV — >• 0. To see this we expand ( !26l) and ( !27|) in fV, 

p = p^ + p. + 0{fV), p=P^0±-{p^ + p_)+O{fV). (28) 
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Without the terms of order fV, this is a solution of the Euler equation ffTSj) 
for constant density p+ + p^. Now the hmit p+ p_ and fV has 
to be taken in such a way that p stays regular and non-negative (we skip 
mathematical subtleties) . 

We define r = {p^/p_)G which, in contrast to the Newtonian case, is 
now restricted to be less than 1. We draw the PSEOS in Fig. ([2]) in terms 
of the rescaled variables p/p-, p/p- and for the same values of r as chosen 
for the NEOS in Fig. namely r = (3 - V5)/2 ^ 0.382, r = 0.6 and 
r = [(v^- l)/2]i/2 ~ 0.786. 



P/P- 




0.2 0.4 0.6 0.8 1 1.2 

P/P- 

Figure 2: The Pant-Sah equation of state ([T]), ([2]) for the values r = 0.382 
(thin solid line) r = 0.6 (medium solid) and r = 0.786 (thick solid line). The 
dotted lines indicate the respective limits of p/p- for p/p- — * oo. 

At first sight the PSEOS looks very different from the Newtonian model. 
Fig. ([T]). In fact, in contrast to the latter, the density now stays bounded 
and tends to p+ + p_ as the pressure goes to infinity (which happens for 
fV 0). This means that for high pressures the PSEOS first violates 
the energy conditions, and finally always becomes infinitely "stiff". Note 
however that Fig. ([1]) and Fig. ([2]) have a very different scale in the p/p- 
direction. For small p/p- and small p/p_ we can still consider the EOS 
as Newtonian limit of the PSEOS (1^ and (ITTI). 

At the surface, the quantity k as defined by ([3]) is related to the speed of 
sound by = dp/dp\s = kJ^- In analogy with ( fT4|) we now determine the 
surface potential, the surface density and Cg from (!25l) . (!26|) and (p71) as 
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fl/-^"^ .-^^ C'-—^ (29) 

^^^-TT7' P'-a+W'' ^^-5(i-r)2- ^2^^ 

Since V^^ < 1, / is bounded from below by / > (1 — r)/(l + r). 

Again ps and are determined by the EOS alone whereas one of / or 
Vs can be used to parametrize the solutions. If r > (6 — a/IT)/5 ~ 0.537, 
then Cs > I, i.e. the speed of sound exceeds the speed of light already at the 
surface. This applies in particular to the curves for r = 0.6 and r = 0.786 of 
Fig. (ED. But r < (6 - v^)/5 implies < 1, and for sufficiently "small" 
spherical stars C < 1 then holds up to the center. An example is the thin line 
T = 0.382 of Fig. ([2]). On the other hand, for all r, C > 1 at the centre if the 
star is sufficiently large due to the stiffness of the PSEOS at high pressures. 
The size limits for the star follow from Fig.l in [TB], and they could also be 
determined from the results of Sect. 3.2. 

While for our purposes the parametric form ( 126|) and ( 1271) suffices as EOS, 
the latter can also be displayed in closed form. We ffist consider the BEOS 
(p+ = 0) which reads 

6/5 

P = 77? ' (30) 

6(py'-pi/5) 

and holds for p < p^. In the general case p+ > it is clearly simplest 
to eliminate one of p_ or p+, and to interpret the other one, together with 
A = fV as parameters of the PSEOS. However, in view of the geometric 
interpretation of p_ and p+, and in view of the "symmetric form" of equations 
I^E^ and (l?n) . it is more natural to eliminate A = fV. To do so we note that 
the following linear combination of equations and ([27]) 

^b(l + A) + 6pAl + ^|,(l-A)-6pAl = 

= ^{l-Xf + ^{l + \r^\' + 2X' + ^ (31) 

gives a polynomial equation of fourth order in A which can be solved alge- 
braically by a standard procedure. Alternatively, we can use ( IHTj) to eliminate 
the fifth and fourth order terms in (I2U1) which leaves us with the polynomial 
equation 

A^ + ^(p + Qp)X^ + ^[1 + 2z^+(p + 5p)]A - z/_p = (32) 

of third order, with 32^^ = p'^^ ± pl^. Solving either (IHTj) or for A = 
A(p, p, p+, pJ) and putting this back again into (^Ul) or (^71) gives the PSEOS 
in closed form H{p, p, p+, p_) = 0. The function H is elementary but involved 
and will not be displayed here. 
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3 The Spherically Symmetric Solutions 



We now determine the spherically symmetric solutions corresponding to the 
two-parameter families of NEOS ([5]) and the PSEOS and ([2]), making use 
of the formulas in the Appendix. By the general theorem |5j there exist 1- 
parameter families of such solutions in either case, for all values of the central 
pressure. We will in particular determine the physically relevant parameters 
mass M and radius R. As well known and easy to see from the definitions and 
the field equations, families of static, spherically symmetric fluid solutions 
are always invariant under the scaling 

p^7V P-^I^P, M^-f'^M, R^'j-^R. (33) 

for any 7 > 0. For our families of solutions this means that one of the 
three parameters is "trivial" in this sense. In Sects. 3.1.2 and 3.2.3 we will 
therefore use scale invariant variables M and R defined by 



M=/-^M, R=f-^R, (34) 

where ps is the surface density. Note that the latter is given in terms of p_ 
and p+ in Newtonian theory by 0141) but in Relativity by 0291) . 

3.1 The Newtonian Solutions 
3.1.1 The matching 

Using Lemma A. 3. of the Appendix with pfj flat, ^ = {v — v)/2, and 
71- = 51271 p_ we can write the spherically symmetric solutions (IMl) of (|TT|) 
as 

It remains to eliminate one of the constants v and /i by global conditions. In 
the case p+ = the NEOS becomes the polytrope of index 5. It follows from 
06 7p that V = which means that JF extends to infinity and 0351) is valid 
for all r. The solutions can be conveniently parametrized by their mass M 
defined in ([9]). 

In the case p+ > 0, it is simplest to parametrize the solutions in terms of 
V which is related to the surface potential by OT^ . To get p it suffices to use 
that f G which implies that 

— - — - -— (36) 
dr dr R 
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where R = r\s is the radius of the star. 

Using and ffT^ and recalhng that v was assumed non-negative, it 
follows that /i can be expressed as 



/i 



and we can write ( l35l) as 



V = < 




M 



TV 



for p+ = 0, 



for p+ > 0. 



(37) 



(3^ 



Note that M can take any value M G [0, oo), and the allowed values for 
the other parameters are v e [0, r] or Vg € [— r, 0]. 

For all p+ and p_ the density, the pressure and the speed of sound follow 
from ffT^ and ffT^ : they are monotonic functions of r. For p+ > the central 
density pc, the central pressure Pc and the speed of sound at the center Cc 
take the values 



Pc 



Pc 



6 



(39) 



These quantities diverge as the parameter v goes to zero. 

Instead of the coordinate expressions (!35l) - (!38l) the matching and the 
solutions can be described in a "covariant" manner in terms oi w = VjfV^f 
which is a function of v in the spherically symmetry case. In particular we 
have w = M~'^v^ in the vacuum region. To determine w for the spheri- 
cally symmetric solutions characterized by TZ_ = 0, we use Lemma A. 2 of 
the Appendix which shows that g~j = gij{v — vY /IQ are spaces of constant 
curvature. With the general formula ( l90l) this yields 



(40) 



Contracting this equation with V*f V-'f and using ([6]) and ( |T2|) gives 



dv 



w 



»7r , 



(41) 
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This has the solution 



w = —p^{v — v)^ [cr^ — [v — t;)^] 



(42) 



for some constant a which has to be determined by global conditions. 

From the exterior form w = M~'^v^ and from the matching conditions ([8]) 
we obtain 

rim rioi} Aoi^ 

(43) 

for p+ > we find 



dw 






dw 




4w, 


dv 


Snp 




dv 


av<s= 


Vs 



Using the asymptotic property ^ for p+ = and 



that 



167rp_M2 



for p+ = 0, 



— for p+ > 0. 

V 



(44) 



Alternatively, equ. (H^ can of course be checked directly from (IHHl) . In 
particular the value of a follows from (1371) or vice versa. 



3.1.2 The Mass-Radius relation 

To determine mass and radius we take equ. (1381) or the first of (H3i) and (jl2 
at the surface and use Vg = —M/R. 



In terms of the rescaled variables 
R^ = -{t-v), 

T 

and implies the mass-radius relation 

M2 - 



this gives 
2v 



(r - vY 



R' 



R.M+ — = 

r 2 



which can be solved for the mass 



M 



tR 



1 ± 



2i?2 



r 



(45) 



(46) 



(47) 



We remark that in ( H71) r could be removed completely by a further 
rescaling of M and R. We avoid this, however, to keep the close anlogy to the 
relativistic case where this is not possible. The behaviour of the parameters 
introduced above is illustrated in Table and Fig. ([3]). 
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Table 1: The parameters of the Newtonian solutions 





V 




R 


M 








dust particle 


T 



















biggest star 


T 

2 


T 

2 


pf 
V 2 


2 V 2 


T(l-2^/2) 
2 




6 


heaviest star 


T 

4 


3r 

4 


1 /i^ 

2V 2 


8 V 2 


7r 

4 


32p_r5 


21p+ 
2 


singularity 





— T 








— oo 


oo 


OO 



0.35 I \ ^ ^ r 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 

R 

Figure 3: The mass-radius diagram for the Newtonian model with equation 
of state ([5]), for the values r = 0.382 (thin line), r = 0.6 (medium) and 
r = 0.786 (thick line). 

For the NEOS ^ with p_ G [0, oo), p+ G (0, oo) the quantities w, v^, Vc, 
Pc and pc can take all values in the open interval bounded by the respective 
values of the "dust particle" and the point singularity, which are clearly 
unphysical themselves. These parameters are monotonic functions of one 
another, and any of them can be used to characterize the solutions. 

On the other hand, the mass and the radius have extrema which follow 
easily from ( H71) and are also given in Table ([T]). Fig. ([3]) where we have 
chosen the same values of r as in Fig. ([T]) shows the following. Starting 
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with the dust particle at R = M = Vs = and increasing pc and Vs, we 
follow the lower branch of the mass-radius curve which corresponds to the 
minus sign in ( H71) . After passing the maximum radius R = \/ the mass 
which is now given by the plus sign in (HTj), continues to increase till the 
"heaviest star" of the table is reached. Then mass and radius drop towards 
the point singularity i? = M = and Vs = —t. The surface potential 
Vs = —M/R (minus the slope of the line joining points of the curve with 
the origin) decreases monotonically along the curve, whence the latter forms 
precisely one "loop". 



3.2 The Relativist ic Solutions 
3.2.1 The matching 

Using Lemma A. 3 of the Appendix with pij = gij, 

and 7l± = 5127rp-|- we write the spherically symmetric solutions of ([2SD as 



1-/^ _ , M±lt^ US) 

g,dx^dx^ = ^^(dr' + r^du^) (49) 

(1 + + 

Again we have to eliminate one of the constants / and fi by global 
conditions. Recall that the parameters in the EOS are now restricted by 
< p+ < p_ < 00 and so < r < 1. In the case p+ = the solutions 
extend to infinity (which can be shown independently of spherical symmetry, 
c.f. [3 EH [25]) and we set / = 1. For p+ > the solutions are finite since 
Ps > 0. We claim that the Buchdahl solutions and the Pant-Sah solutions 
are given by (HH!) and (H9|) with 



for p+ = 0, 



IJ= { ^M\jTip_ (50) 
for p+ > (J, 



and 



r\lil±ry+ ^ ^ ^ (51) 



which requires / < 1 to make sense. For p_|_ = equ. (1501) follows easily from 
the asymptotic condition (l20l) . On the other hand, for p_|_ > the matching 
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to Schwarzschild is quite involved as the isotropic coordinates of 
which simphfy the interior solutions are unsuited for the matching. We will 
verify ( !50l) below by matching " covariantly" (c.f. ( !59l) ). 

We also write the Pant-Sah solutions in the alternative form 



where 



ds^ = -V'dt^ + n\dr^ + r^cLj') 



2/x2 - S2 
and the constant S is given by 



[i + fvy ii-fv)\ 



/i for = 0, 
for p+ > 



(52) 
(53) 



(54) 



in terms of the quantities defined in fl50p and (1511) . The form (1521) will be 
useful in Section 5 for the proof of spherical symmetry of solutions with the 
PSEOS. Equ. ( |53l) makes sense for p+ = as well (i.e. for the BEOS with 
p_ 7^ as well as for vacuum p_ = 0) and reduces to = 4/(1 + V)^ in 
either case. 

In analogy with the quantity w in the Newtonian case, we now determine 
W = V^WiV which is a function of V in spherical symmetry. In particular, 
for Schwarzschild we have W = (16M)~2 (i _ V"2) . For our model charac- 
terized by TZ± = const., Lemma A. 2 shows that the spherically symmetric 
solutions are spaces of constant curvature. Using the field equations (|T71) and 
the general formula (pUj) with pij = gij, = (1 ± fV)/2, pij = gf- we have 

= V(l ± VfBf. = C[(l - fV^)ViVjV + 6fVViVV,V]. (55) 



Contracting this equation with and using the field equation (|T6 

and (EED and (EZD gives 



_d_ 
dV 



W 



_(l-/2\/2)4 

which has the solution 
Tip 



P- 



47r(l - fV) 

37(1+77) [{i + fvy {i~fvy\ 



w 



3/2 

7rp_(2p2 



1 - fvY 

S2^ 



[i + fvy 



12/2 



1 - fvyn 



;i + fv) 



(56) 

(57) 
(58) 
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In equ. flFTI) T? arises as a constant of integration, and we first verify that 
it is consistent witfi tfie earlier definitions (jSll), flHIj) and fj^ . 

For p+ = this follows once again from the asymptotics, equ. (l20l) . For 
p+ > this is done with the "covariant" matching condition ( fT9l) which 
becomes 



dV 





'dW 








_dV _ 







(59) 



Next, with a little algebra one can check (!3H|) which contains and f2 defined 
in flSOp and (1531) . To verify that /i as defined in (1501) in fact agrees with the 
constant appearing in (l48l) it is simplest to use (|57j) and the general definition 
W = ViVV'V. 

Finally we note that one can alternatively write the Pant-Sah solutions 
by using \^ as a coordinate everywhere. (Equ. (13^ still contains "r"). From 
eqs. (3.1) and (3.17) of [12] one finds [7j 



^47rV-(2/i2-S2)2(l-/2K2)6 



du' 



(60) 



3.2.2 The centre 

We now turn to the important issue of regularity at the centre. The latter 
is characterized either by r = or by the minimum of V, i.e. W = 0. From 
(HHl) and (H9l) it is easy to see that the centre is regular if > 0; it can 
be made manifestly regular, i.e. Qijdx^dx^ = (1 + 0{r'^)){dr'^ + r'^duj'^) by a 
suitable rescaling of r. 

We also note that either from fj^Hl) or fl3Sl) it follows that regularity is 
equivalent to /i < 1. For p+ = this entails the lower hound > 3/167rp_ 
/or t/ie mass, while there is no upper bound. For p_|_ > we have collected 
in Table ([2]) the most important parameters which are monotonic functions 
of one another and provide unique characterizations of the model. We use 
the shorthand T± = Vl ± r + r^. 

Table 2: Some parameters for the Pant-Sah solutions 





/ 


Vs 




pc 


Pc 


dust particle 


1 
1 


— T 

+r 


1 


1 


l+r4 





singular centre 


(1- 
(1+ 









P- + P+ 


oo 
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The allowed parameter values are bounded by the respective ones of the 
"dust particle" with pc = and Vc = Vg = 1, and the model with singular 
centre for which pc = oo and Vc = 0. Like their Newtonian counterparts these 
limits are unphysical, but unlike the Newtonian ones the singular solution 
has now finite extent. The dust particle has pc = Ps > 0, while as the singular 
centre is approached pc always stays finite (in contrast to the Newtonian case) 
due to the "stiffness" of the PSEOS. This singular model also has largest 
redshift, which can be tested against the Buchdahl limit [33] Vg > 1/3. 
The latter is saturated for fluids of constant density only. Such fluids are 
approached by the present models for r ^ 1 (c.f. Equ.(|28l)). and in fact we 
flnd that K — 1/3 in this limit. 



3.2.3 The mass-radius relation 

To obtain mass and radius we use 1 — = 2M/R, fl3Ul) and flFTI) at the 
surface. In terms of the rescaled variables (IMl) we obtain 



2 r 



1-/ 
At 

16r(l + r)4 



P 



:i+r) 



Eliminating / gives the mass-radius relation 



2r 



P 



R' 



[1 + R')RM + — = 



which can be solved for the mass 



(61) 
(62) 



(63) 



M 



tR 



[l + r) 



1 + R^± 



(64) 



The extrema of mass and radius are listed in Table ([3]). (Recall that 
T± = Vl±r + r2). 

As in the Newtonian case the surface potential characterizes the solution 
uniquely, which implies the loop-like structure of the mass-radius curves, Fig. 
(jlj). We flrst describe the diagrams for sufficiently small values of r such as 
r = (3 — V5)/2 ^ 0.382 and r = 0.6. (These particular values correspond 
to Vg = 1/3 and V^^ = 1/6 at the respective mass maxima). Starting with 
the dust particle and increasing p^, K decreases and the mass-radius curve 
corresponds to the minus sign in front of the root in flM|) . At the maximum 
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Table 3: Surface potential, radius and mass of the Pant-Sali solution 





V"2 

s 


R 


M 


dust particle 


1 








biggest star 


1-T 
1+T 




1+T 


heaviest star 


(1-T)(2r+ + 1-T) 


r++T-i 


2Tf-(2+T)(l+2T)(l-T) 


3(1+t)2 




3V3^(1+t)2 


sing, centre 




2tV^(t2+1) 


4T2(l+r2)^r(T2+l) 






r2rf 



0.4 




0.35 




0.3 




0.25 




M 0.2 




0.15 




0.1 




0.05 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

R 

Figure 4: The mass-radius diagram for the Pant-Sah solution for the values 
r = 0.382 (thin line), r = 0.6 (medium) and r = 0.786 (thick line). 

radius which is now at i? = i/r, we pass to the plus sign. From then onwards 
the star shrinks, reaching its maxium mass at some lower value of Vg, and 
subsequently losing weight as well. In contrast to the Newtonian case, the 
singular model now prevents the "mass-radius loop" from closing. As already 
mentioned in the discussion of Table ([2]), at some finite size of the star the 
central pressure diverges, and this is where the curves in Fig. (jlj) terminate. 

For T = 0.6, the star with maximal mass still has a regular centre. How- 
ever, for larger values of r, the central pressure diverges before the mass 
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maximum or even the maximal radius are reached. This means that the 
"biggest star" and the "heaviest star" in Table (3) only make sense if the re- 
spective values of Vs are larger than the ones given for the " singular centre" . 
For r = [(-\/5 — l)/2]2 ^ 0.786 the star with maximal radius is precisely the 
first one with singular centre, and the meaningful part of the mass-radius 
curve is monotonic (c.f. Fig. (jlj)). 

We finally note that the mass radius curves for the "softer" PSEOS (such 
as r = 0.382 in Fig (jlj)) resemble strikingly the mass-radius curves for quark 
stars [19]-[22] (in particular those with "harder" equations of state such as 
[21], [22]). Putting p„ = 3 GeV/fm^ in the PSEOS with r = 0.382 one 
obtains typical values of about 1.5 solar masses for the maximal mass and 
about 7 km for the maximal radius. However, this coincidence should not 
be overestimated. As mass and radius are obtained by integration, they 
"smooth out" differences in the EOS, which seem quite substantial even at 
moderate densities. Moreover, we recall that for the PSEOS the pressure al- 
ways diverges at finite density, whereas the EOS of "ultrarelativistic" quarks 
is nowhere too far from p = p/3. This discrepancy prevents us from mod- 
elling extreme quark conditions and has rather drastic consequences for the 
mass-radius relation. As follows from Harrison et. al. [T7] and has been 
shown rigorously by Makino [18], if the quotient p{p)/p for some given EOS 
tends to a constant sufficiently fast for p —>■ oo or p ^ oo, the mass-radius 
curve develops the form of a " spiral" , with an infinite number of twists, for 
high central pressure. While the EOS for quark stars given in the literature 
seem safely within the range of the Makino theorem, the mass-radius dia- 
grams are normally not drawn till the spiral sets on. On the other hand, the 
mass-radius diagrams for the Pant-Sah solutions are single open loops, which 
we have drawn to the end (infinite central pressure) in Fig. (jl]). 

4 Proofs of Spherical Symmetry 

In Newtonian theory Lichtenstein [35] has given a proof of spherical sym- 
metry of static perfect fluids which satisfy p > and p > 0. Under the 
same condition on the equation of state, Masood-ul-Alam has recently proven 
spherical symmetry in the relativistic case by using a substantial extension 
of the positive mass theorem [3j. For the relativistic model considered in 
this paper, spherical symmetry is a consequence of the uniqueness theorem 
of Beig and Simon [12]. In Sect. 4.2 we reproduce the core of this proof for 
the present model, for which it simplifies. 

In Sect. 4.1 we give a version of the Newtonian proof which resembles as 
closely as possible the relativistic proof, substituting the positive mass the- 
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orem by the virial theorem. A proof along the same hnes has been sketched 
in [12] for fluids of constant density. 



4.1 The Newtonian Case 

We use the notation of sections 2.1.1 and 2.1.2 with the following additions 
and modification. We define w = g^^VivVjV as in Sect. 2.1.1. However, for 
a given model described by Vg, we now denote by wo{v) the function of v and 
Vg defined by the r.h. side equ. (H2l) . Note that this function may become 
negative, which happens if the central potential Vc of the given model is less 
than the central potential of the spherical symmetric model with the same 
Vg- The proof of spherical symmetry proceeds by showing that w and Wq 
coincide. We split this demonstation into a series of Lemmas. 



Lemma 4.1.1. (The virial theorem) 

For all static asymptotically flat fluids as described in Sect. 2.1.1., we have, 
denoting the volume element by drj, 

(6p + pv)drj = 0. (65) 

Remark. In kinetic gas theory, the two terms in the integral flU31) are four 
times the kinetic energy and twice the potential energy of the particles, re- 
spectively. 
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Proof. Let C,i be a dilation in fiat space, i.e. V(j^j) = gij. (In cartesian 
coordinates, gij = 6ij and = x*). Let v, p and p define an asymptotically 
flat Newtonian model as in Sect. 2.1.1.. Then there holds the Pohozaev 
identity [37] 



v. 



e V.i; + -V 



27r(6p + pv) 



(66) 



which is verified easily. We integrate this equation over Ai and apply the 
divergence theorem, using that the integrand in brackets on the left is contin- 
uous at the surface. Due to the asymptotic conditions ©, ffTUj) the boundary 
integral at infinity vanishes which gives the required result flU31) . □ 



Lemma 4.1.2. For solutions with the NEOS ([5]), we have 
= j {6p + pv)dr] = -p+Y + vM 



= iZ^d, 



{v - vY 



(67) 
(68) 
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where Y = jjrdv is the volume of the fluid. In particular, p_|_ = iff the 
solutions extend to infinity {0 = Vg = v — r) 

Proof. For the Newtonian model the virial theorem ( l65i) and ( fT2l) imply 

= y {6p + pv)dr] = J [6p — p{v — v)]dr] + V J pdr] = —p^Y + vM (69) 

which proves (!U7|) . To show (lUHl) we use the divergence theorem, IQ, (fT^ . 
([HD and (112D. We obtain 



W-Wq 



dr] 



= 1/ v..d5'-i=^^/*, = l^fM-^!±I) (70) 

Jd:F V 'T-^ \ V J 

Equ. (1U7|1 and (ITUjl now give the required result. □ 
Lemma 4.1.3. For the NEOS there holds, inside T 

Proof. We first note that formula fpl|) with = (7^^, $ = (-u — f and 
pij = ^^pij = Qij implies B~, = —{v — v)^^C-[V^Vj {v — v)"^]. If follows that 



2A-J-^ = -Vl[B-^VUv - vf] = [v - v)%B'l - 

- ^[VUt; - vr]Vrn- = iv- vfBr^B'l (72) 

where we have used the Bianchi identity V_Blj = VjTZ^ /6 and the fact that 
7^- = const, for the NEOS. □ 

Lemma 4.1.4. In V we have 

\v\'^[v — vj [v — 

where Bij is the trace-free part of the Ricci tensor w.r. to the metric 
9ij = v^gij- 
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Proof. In the vacuum region fl7T]) still holds (we set f = 0), since the 
metric 'gij has curvature 7^ = 0. It follows that 



A- 



W — Wq 



V 



\v\-^{v — v) {v — v)^ 
{v — v)^ 



A 



[V — V) 
W — Wq 



\v\ 



W — Wq 



{v — v) 



V — V 
7 



where V* and A refer to ^ij. 

Lemma 4.1.5. On = V U JF, we obtain w < wq. 



(74) 



□ 



Proof. The weak maximum principle applied to f[7T|) on JF implies that 
{w — Wq) /{v — v)^ takes on its maximum at some point q G dT, i.e. 



W — Wq w 

sup — — < max ■ 



Wq 



W — Wq 



[V — V] 



dT (V — V) 



(75) 



On the other hand, the weak maximum principle applied to ( 1731) shows that 
either {w — wq)/\v\^{v — v) takes on its (absolute) maximum at infinity (where 
it vanishes) or that it has a positive maximum on dJ-'. In the latter case 
this maximum is located at q since v is constant on dJ-". This leads to a 
contradiction as follows. Taking n* to be the normal to dJ^ directed towards 
infinity, we have 



W — Wq 



{v — vY 



W — Wq 



{v — v)\v\^ 



+ 3 



V[W — Wq) 



n^V. I^IL. (76) 



By the boundary point lemma |36j applied to (I7H1) . the first term on the 
right of flTHj) is negative, and the same applies to the second term by virtue 
of Av = 0. It follows that 



W — Wq 



V 



V 



< 



(77) 



But as V — V is on Ai, this contradicts ( 1751) . Hence we are left with 
{w — WQ)/\vf{v — v) < in V which, together with ( 1751) implies {w — wq)/{v — 
f )^ < everywhere on A4. This proves the Lemma. □ 

Lemma 4.1.6. On Ai = VUJ-',we have w = Wq. Furthermore w = Wq{v) is 
positive for v > Vmin, smooth in [vmin, Vs) and such that at Vmin there holds 
Wq = and dw/dv ^ 8np. 
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Proof. The first assertion follows from Lemmas 4.1.2 and 4.1.5., and tlie 



Proposition 4.1.7. Asymptotically fiat solutions with the NEOS are spher- 
ically symmetric and uniquely defined by Wq. 

Proof. A trivial modification of a relativistic result, Lemma 4 of [12], has the 
conclusion of Lemma 4.1.6. as hypothesis. The conclusion of this modified 
Lemma is the proposition. 

4.2 The Relativistic Case 

We use the notation of the previous sections with modifications analogous 
to the Newtonian case. We recall that W = g^^ViW jV, and for a given 
model described by Vg we now denote by lyo(^) the function of V and Vg 
defined by the r.h. side equ. flFTj) . Again this function becomes negative if 
the central potential Vc of the given model is less than the central potential 
of the spherical symmetric model with the same Vg. We first prove that W 
and Wo coincide [27], which is done in a series of Lemmas. From Lemma 
4.2.3. onwards, they are direct counterparts of the Newtonian ones in the 
previous section. 

Lemma 4.2.1. (The vanishing mass theorem). We recall that an 
asymptotically fiat Riemannian manifold with non-negative scalar curvature 
and vanishing mass is flat |26j . 

Lemma 4.2.2. For fluids with PSEOS, we have 



defined in f l3Hl) . For p+ = 0, i.e. for the BEOS and for vacuum, we obtain 

n, = o. 

Proof. For the curvature w.r.t. to g*^ we obtain 



rest is easily checked. 



□ 





) 



2 



7^. 



2 



2n 



dV^ 



{Wo - W) 



167rfi2 p + (p + 3p) 



V dVt 

n~dv 



( 



1 



W_ 

Wo 



) 



(79) 
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Here the first equation holds for conformal rescalings of the form g*j = Q''^gij 
(for any Qij and ^(V^) if 7?. is a function of V only) while the second one uses 
the general formula ( !89|) , property ( l52l) , and the field equations (fT6l) and ( fT71) . 
Now (ITHI) follows by using the explicit forms ( l26l) . ( 1271) and ( l53i) of p, p and fl. 

Lemma 4.2.3. For solutions with PSEOS, we have 

W-Wo _ V\l±fVl ^ , 
^ (1 _ /2v^2)4 - (1 ^ /^)io ^^.^± ^ (80) 

where A' refers to the metric g'^j = (1 — f'^V^)^gij/16V'^. 

Proof. We first note that formula ^ with pij = gf^, $± = 2/(1 ± fV), 
and pij = ^±pij = gij implies, together with the field equation ffTTj) . 



[V±Xf ] = a^Bt, (81) 

where we have defined 
Then we find from 08 ip that 

— — ^ (1 _ /2y2)4 - ^± L%^±J - « + 

+ ^XiVfn^ = a^Bpt (83) 
where we have used the Bianchi identity ^±Bf, = and the fact that 

I J J 

TZ = const, for our model. □ 

Note that the argument of the Laplacian on the l.h. side of (IHOl) agrees 
with TZ-i, as given in (!7S|) modulo a constant factor. In other words, (!75|) and 
show that A'7^, > 0. 



Lemma 4.2.4. In V, we have 



where Bf^^ are the trace free parts of the Ricci tensors w.r. to the metrics 



^5 = (1 ± V)%,/2. 
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Proof. In vacuum flHUl) still holds (we set / = 1), since the metrics ^ 
have vanishing curvatures TZ^. It follows that 

W-Wo 



A' 



[1 - v^y^i - pv^) 



2\6 



2t/2\2 



± 



(l-\/2)2 
W -Wo 



1 _ y2 



1-/2^2 (l_y2)4 



^5) 
□ 



Proof. The weak maximum principle applied to flHOj) on JF implies that 
(ly — iyo)/(l ~ /2\/2^4 ^g^jjgg its maximum on some point q G i.e. 



(l-/2\/2)5" (1-1/2)4 (l-/2V2)5(l^y)5^ij^±- 

where and refer to respectively. 
Lemma 4.2.5. On M = V U dJ^ U we have W < Wq. 



W -Wq W -Wq _ W -W i 

'7 (1 _ /2y2)4 - (1 _ /2V^2)4 - (1 _ /21/ 







2T/2^4 



(86) 



On the other hand, the weak maximum principle applied to ( l84l) shows that 
either {W — Wo)/{l — V'^)^{1 — pV"^) takes on its (absolute) maximum at 
infinity (where it vanishes) or that it has a positive maximum at some point 
q G ST. In the latter case this maximum is located at g, as V is constant on 
(9JF. This leads to a contradicition as follows. Taking n* to be the normal to 
(9JF directed towards infinity, we have 



W -Wo 



(1 - pv^y 



(l + r)2(l-KT 



W -Wo 



{I- PV^){l-V^) 



?,{l + Tf{l-p)V,{W-Wo 
32r3(l - K,2)3 



n'Vi V\ 



17) 



By the boundary point lemma applied to fl84p . the first term on the right 
of 0871) is negative, while the second term (without the minus in front) is 
positive by virtue of AV = 0. It follows that 



W -Wo 



pv'^y 



< 



But as 1 — /2V"2 is on A^, this contradicts ( |88l) . Hence we are left with 
{W - Wo)/ll - ^2)3(1 - /2-^2) < in V which, together with implies 
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[W — Wo)/{l — pV'^Y — everywhere on Ai. This proves the Lemma. □ 

Lemma 4.2.6. On = V U JF we have W = Wq. Furthermore W = Wo{v) 
is positive for V > Vmin, smooth in [Vmin, and such that at Vmm there 
holds Wo = and dW/dV ^ 87r(p + 3p)V. 

Proof. The first assertion follows from Lemmas 4.2.2 and 4.2.5., and the 
rest is easily checked. □. 

Proposition 4.2.7. Asymptotically flat solutions with the PSEOS are 
spherically symmetric and uniquely defiend by Wq. 

Proof. Lemma 4 of [12] (which is a reformulation of results of [3H1 EH]) 
has the conclusion of Lemma 4.2.6 as hypothesis. The conclusion of the for- 
mer Lemma is the proposition. 

Acknowledgements. I am grateful to J. Mark Heinzle, Kayll Lake, 
Marc Mars and to the referees for useful comments on the manuscript. 

5 Appendix 

We recall here the basic formulas for the behaviour of the curvature under 
conformal rescalings of the metric and the standard form of metrics of con- 
stant curvature. We give the proofs of the latter two Lemmas. 

Lemma A.l. Let pij, pij = ^'^pij be conformally related metrics on a 
3-dimensional manifold A4. Then the scalar curvatures 3?, 9fJ and the trace- 
free parts Bij = C[^ij] and Bij = C[^ij] of the Ricci tensors 3ftjj and 
behave as 



where Vj is the gradient and A = p'-' VjVj the Laplacian of pij. 

Lemma A. 2. Any smooth, spherically symmetric metric {A4,pij) with 
constant scalar curvature 3? is a space of constant curvature (i.e. Bij = 0) 
and can be written as 




(89) 




B^j - 2<!?-^C[V^Vj^ + 6$"^C[Vi$Vj$] 



(90) 



ds^ 



1 



^{dr^ + Pdu^) 



(91) 



(1 + ^5?r2) 
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Proof. By solving ODEs we can write pij in isotropic coordinates as 

ds^ = ^{r)\dr^ + r^du^) (92) 

To determine $, we solve (1891) with pij flat and 3? = const.. Near the 
center, the solution determined uniquely by the initial values <l'c = 1 and 
d^/dx^c = reads 

$(r) = ^ (93) 

^1 + ^3ftr2 



24" 

The solution can be extended analytically as far as required, which gives 

dnn). □ 

Lemma A. 3. Let (A^+, p^) and (A^_,p^) be 3- dimensional, spherically 
symmetric manifolds with smooth metrics and with constant scalar curva- 
tures 3fJ+,3fJ_. Then the smooth, spherically symmetric solutions $+(r) to 
the equation flHUj) on (A^+, pfj) are given by 



1 + J-gf? , 



where /i is a constant. 



Proof. Using Lemma A. 2., the required solution is determined by a con- 
formal rescaling p^- = between spaces of constant curvature. Writing 
the metrics in the forms fpij) gives 

'^+^''+^' J drl + rlduj^) = .{dr^ + duo^) (95) 

Hence ^+(r+) defined by 




with a yet unknown relation r_ = r_(r+), determines a conformal diffeomor- 
phism of flat space, i.e. 

^X{drl + rlduo^) = {drl + r^doo^) (97) 

By (l89l) all such diffeomorphisms are solutions of A\E'+ = on flat space, 
and hence given by \E'+ = fi + u/r^ for some constants fi and z/. A consistency 
check with f l97p now shows that either u = and r_ = /i^r+ or = and 
r_ = z/^/r+. Setting r+ = r, the flrst case leads directly to while in 
the second case we have to put /i^ = 24/(z/^$ft_). The solution again extends 
analytically as far as needed. □ 
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